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ABSTRACT 

Moire  interferometry  was  used  to  determine  the  two  orthogonal  displacement  fields 
surrounding  a  stably  growing  crack  in  a  thin  2024-T3  aluminum  alloy,  single-edge  notched 
(SEN)  specimen.  The  measured  displacements  were  used  to  compute  the  J-integral  and 
the  associated  HRR  and  the  second  order  Q  displacement  fields.  Also  a  2-D  elastic-plastic 
finite  element  model  of  the  fracturing  specimen  was  executed  in  its  generation  mode  to 
compute  the  crack  tip  displacement  field  and  the  J-integral  for  small  crack  extension.  The 
numerically  computed  and  experimentally  derived  J-values  agreed  prior  to  crack 
extension  but  differed  increasingly  with  stable  crack  growth.  The  HRR  and  measured  v- 
displacements  agreed  well  but  significant  difference  between  the  HRR  and  measured  u- 
displacements  was  noted.  Thus  the  Q  component  of  the  v-displacement  was  negligible  but 
the  corresponding  u-component  varied  nonlinearly  with  radial  distance  from  the  crack  tip. 

INTRODUCTION 

A  hybrid  experimental-numerical  procedure  using  the  displacement  field  obtained  by 
geometric  moire  technique  to  drive  elastic-plastic,  2-D  and  3-D  finite  element  codes  was 
first  demonstrated  by  Hareesh  and  Chiang  [1].  In  their  analysis,  experimental  surface 
displacements  were  used  to  compute  the  interior  two-  and  three-dimensional 
displacements  and  strains  in  a  plastically  deforming  stationary  notch  tip  of  a  3.2-mm  thick, 
6061-T6  aluminum  alloy,  single  edge  notched  (SEN)  specimen.  The  authors  concluded  that 
a  ten  fold  saving  in  computer  processing  time  was  achieved  through  their  hybrid  analysis 
and  that  their  3-D  hybrid  model  showed  significant  three  dimensional  effects,  which 
otherwise  could  not  be  detected  by  surface  displacement  measurements  alone,  in  the  very 
vicinity  of  the  notch  tip. 

With  the  increasing  use  of  an  order-of-magnitude  more  sensitive  moire 
interferometry,  it  was  inevitable  that  it  be  incorporated  into  the  above  hybrid  analysis. 
Sivaneri  et  al  [2]  used  the  displacement  fields  obtained  by  a  white-light  moire 
interferometry  of  a  2024-0,  0.8-mm  thick,  fatigue  precracked  SEN  specimen  [3]  to  drive  an 
elastic-plastic,  plane  stress  finite  element  code.  Among  other  numerical  results,  the  J- 
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integral  value  and  the  associated  HRR  displacement  field  [4,5]  at  the  initial  stage  of 
loading  were  computed.  The  latter  was  in  excellent  agreement  with  the  two  orthogonal 
measured  displacements. 

In  recent  years,  one  of  the  authors  and  his  colleagues  have  been  using  the 
orthogonal,  crack  tip  displacement  field  provided  by  moire  interferometry  to  determine 
the  J-integral  value  and  the  associated  HRR  crack  tip  displacement  field  surrounding  a 
stably  growing  crack  in  thin,  2024-0,  2024-T3,  5052-H32  and  2091 -T3  aluminum  alloy  SEN 
specimens  (6-8).  These  experimental  results  showed  that  these  J-integral  values  differed 
substantially  with  previously  published  solutions  [9,10]  in  the  presence  of  large  scale 
yielding  and  were  not  path  independent  after  a  moderate  stable  crack  growth  of  5  mm. 
Also  the  HRR  displacement  field,  which  was  computed  by  using  the  measured  J-integral 
value,  did  not  agree  with  the  measured  displacement  field  for  submillimeter  crack  growth. 
These  conclusions  were  unexpected  since  both  the  J-integral  and  the  HRR  field  had  been 
studied  in  detail  and  had  been  repeatedly  verified  through  many  numerical  investigations. 
These  results  also  suggest  that  perhaps  the  fault  lies  in  the  elastic-plastic  finite  element 
codes  based  on  deformation /incremental  plasticity  with  isotropic/kinematic  strain 
hardening.  Furthermore,  the  slope  of  the  log-log  plots  of  the  second  order  Q  displacement 
components  varied  irregularly  with  no  consistent  second  order  strain  singularity. 

The  purpose  of  this  paper  is  to  compare  the  elastic-plastic  fracture  parameters, 
obtained  from  a  commercially  available  2D  finite  element  code  with  those  obtained 
experimentally  surrounding  a  stably  growing  crack  in  a  thin  fracture  specimen.  Since  the 
existing  2D  finite  element  codes  generally  predict  the  overall  elastic-plastic  response  of  a 
structure  with  reasonable  accuracy,  this  study  focuses  on  the  ability  of  a  code  to  accurately 
compute  the  J-integral  and  the  crack  tip  displacements.  Thus  the  FEM  code  will  be 
confined  to  a  local  region,  surrounding  the  crack  tip,  by  prescribing  the  measured 
displacements  along  a  remote  boundary  of  this  stationary  region. 


TWO-PARAMETER  CRACK  TIP  STRESS  HELD 

The  two-parameter  representation  of  the  elastic-plastic  crack  tip  stress  field  has 
been  discussed  by  Li  and  Wang  [11],  Sharma  and  Aravas  [12]  and  Yang,  Cho  and  Sutton 
[13].  This  plane  strain,  asymptotic  stress  field  is  based  on  the  J2  deformation  theory  with 
the  following  power  hardening  material  of  Ramberg-Osgood: 
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where  i,  j  =  1  or  2  corresponds  to  a  Cartesian  coordinate  system  with  axes  parallel  or 
perpendicular  to  the  crack,  respectively,  E  and  n  are  the  modulus  of  elasticity  and  Poisson's 
ratio  respectively,  Oo  and  e©  are  the  yield  stress  and  strain,  respectively,  a  and  n  are 
material  constants  and  Ceand  Sij  are  the  equivalent  and  deviatoric  stresses,  respertively 
The  existence  of  the  following  asymptotic  expansion  of  the  crack  tip  solution  was  then 
postulated  and  represented  in  terms  of  the  equivalent  stress  as 

=  rso«(o)(e)  +  ft  Oe(l)(0)  +  •  as  r  0  (2) 
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A  J-integral  is  then  used  to  determine  the  leading  order  exponent  of  s  =  -l/(n+l)  and 
to  no  surprise  the  HRR  crack  tip  field  [4]  is  recovered  from  the  first  term  in  Equation  (2). 
The  order  of  magnitude  of  the  second  term  in  Equation  (2)  depends  on  the  second 
exponent,  t.  In  particular,  when  t  =  0  for  the  second  and  all  higher  order  terms,  the 
simplified  Q  stress  component  of  [14,15]  is  recovered. 

Based  on  extensive  numerical  analysis,  O'Dowd  and  Shih  [14]  suggested  the 
following  two-term  approximation  to  Equation  (2)  in  terms  of  a  polar  coordinate  at  the 
crack  tip  as 


J_ 

f  rirr  GrQ  ^ 

f  OrB  '' 

Oo 

t  CTre  009  J 

[^(X£o<Tolnr  J 

I 

dee  J 

^  (j/aoJ  ^ 

I 

Oee  J 


(3) 


where  J  is  Rice's  J-integral,  Q  is  a  dimensionless  constant  that  controls  the  magnitude  of 
the  second  stress  term  and  In  =  In  (q)  is  tabulated  in  Reference  [16].  Equation  (3)  is 
specifically  restricted  to  the  forward  sector  of  the  crack  tip. 


O'Dowd  and  Shih  further  reasoned  that  since  dn-  =  Oee  and  I  Ore  I  «  I  ^ee  I ,  Q  is 
essentially  a  plane  strain,  stress  triaxiality  parameter  which  could  include  the  influence  of 
all  higher  order  terms  in  Equation  (2)  [14].  The  strain  and  displacement  components 
corresponding  to  this  plane  strain  Q  component  are  then  a  hydrostatic  strain  and  a 
linearly  varying  homogeneous  displacement,  respectively. 

For  the  state  of  plane  stress.  Reference  [12]  shows  that  the  second  order  stress 
component  of  in  Equation  (2)  approaches  infinity  as  0  — >  160°  and  thus  the  second  order 

solution  probably  is  not  separable.  In  the  region  of  0  <  140°,  however,  Gg  ^  variation  is 

normal  and  thus  one  can  speculate  that  the  two  parameter  stress  expression  holds  within 
this  restriction,  or  specifically  in  the  forward  section  of  the  crack  tip  as  noted  in  Reference 
[14].  This  hypothesis  is  important  if  the  two-parameter  J-Q  theory  is  to  be  used  in 
characterizing  ductile  fracture  of  thin  plates  where  a  100  percent  shear  lip  with  no 
cleavage  fracture  is  anticipated. 

The  strain  and  the  displacement  components  corresponding  to  the  above  two  term 
representation  of  Equation  (3)  can  be  represented  in  a  nondimensional  form  as  [14] 
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where  u^i  is  the  rigid  body  displacement. 

In  the  following,  the  possibility  of  extending  the  plane  strain,  two-parameter 
characterization  of  a  stationary  ductile  crack  to  large  stable  crack  growth  problems  in  thin 
plates  is  investigated.  Specifically,  the  experimentally  determined  crack  tip  displacement 
field  in  a  thin  plate  specimen  with  that  predicted  by  the  J-Q  theory  is  compared.  Only  a 
qualitative  comparison  is  possible  since  the  second  order  term  in  Equation  (6)  could  not  be 

computed  as  u^^^  (q)  is  not  known  for  the  state  of  plane  stress. 

EXPERIMENTAL  ANALYSIS 
Experimental  Procedure 

The  orthogonal  displacement  components  surrounding  the  crack  tip  were  measured 
by  moire  interferometry  using  relatively  coarse,  cross  diffraction  grating  of  40  lines/mm. 
This  coarse  grating  was  necessary  due  to  the  gross  yielding  and  the  large  strain 
components,  which  will  generate  a  moire  fringe  pattern  too  dense  to  resolve,  associated 
with  the  large  stable  crack  growth  in  a  ductile  specimen.  The  coarse,  crossed  diffraction 
grating  was  transferred  onto  the  specimen  surface  using  photoresist  and  is  similar  to  the 
procedure  developed  by  Ifju  and  Post  [17].  However,  in  this  study,  the  highly  polished 
surface  of  the  aluminum  specimen  provided  sufficient  reflectivity  and  thus  an  evaporated 
aluminized  coating  was  not  used.  This  reflective  specimen  surface  also  eliminates  the  loss 
of  moire  fringes  at  high  strain  where  an  aluminized  coating  craze  and  obliterate  the 
diffraction  grating. 

The  specimen  was  then  illuminated  by  a  four  beam  moire  interferometer  [18]  for 
simultaneous  recording  of  the  two  orthogonal  displacement  fields.  Figure  1  shows  the 
moire  interferometry  setup  and  the  u-v  mirror  arrangement  used  in  this  study.  The  coarse 
diffraction  grating  reduced  the  incident  angle  of  the  four  beams  thus  simplifying  the  u-v 
mirror  supports. 


Specimen 

The  specimen  consists  of  a  fatigue  precracked,  thin  single-edged  notch  (SEN),  2024- 
T3  aluminum  alloy  specimen  shown  in  Figure  2.  The  moire  diffraction  grating  covered  a 
region  of  25.4  x  50.8  mm  surrounding  the  crack  as  shown.  The  SEN  specimen  was 
subjected  to  uniaxial  tensile  loading  in  a  displacement  controlled  testing  machine  and  the 
moire  interferometry  patterns  were  recorded  at  various  stages  of  stable  crack  growth. 


NUMERICAL  ANALYSIS 

A  commercial  finite  element  code  was  used  to  compute  the  elastic-plastic  state 
associated  with  stable  crack  growth  in  this  specimen.  The  objective  of  this  numerical 
analysis  was  to  generate  numerical  results,  which  can  be  compared  with  the  experimental 
results,  at  various  stages  of  stable  crack  growth  with  large  scale  yielding.  In  practice 
however,  the  available  code  could  only  handle  modest  plastic  straining  associated  with  a 
small  stable  crack  growth  of  Aa  =  1.5  mm  before  collapsing. 

Unlike  traditional  finite  element  method  (FEM),  the  measured  displacements  near 
the  boundary  of  moire  grating  were  used  as  input  boundary  conditions  to  the  FEM  model 

4 


of  the  SEN  specimen.  As  described  in  [1,2],  this  procedure  not  only  results  in  saved 
computer  time  but  provides  detailed  information,  which  is  lacking  in  the  moire  analysis 
particularly  at  the  initial  stage  of  stable  crack  growth,  in  the  immediate  vicinity  of  the 
crack  tip.  Also,  as  described  in  [1,2],  this  load  path  dependent,  elastic-plastic  finite  element 
analysis  must  start  from  the  early  stage  of  plastic  yielding  and  proceed  incrementally 
along  the  loading  path. 

Figure  3  shows  the  finite  element  model  used  in  this  analysis.  The  power  hardening 
stress-strain  relation  used  in  this  incremental  elastic-plastic  analysis  was  that  of  Dadkhah 
[19].  Also  due  to  the  sensitivity  of  the  FEM  to  a  displacement-prescribed  boundary 
condition,  a  second  order  curve,  which  was  fitted  to  the  measured  boundary  displacements 
obtained  from  moire  interferometry,  was  used  as  the  input  boundary  condition. 

RESULTS 

A  total  of  nine  increments  of  load,  as  shown  in  Figure  4,  were  applied  and  the 
corresponding  moire  interferometry  fringes  were  recorded.  Figure  5  shows  typical  moire 
fringe  patterns  corresponding  to  the  displacement  parallel,  u,  and  perpendicular,  v,  to  the 
crack.  Also  shown  in  Figure  5  are  the  rectangular  contours  used  for  the  J-integral 
computation  [19].  This  computation,  using  moire  interferometry  data,  was  conducted  only 
fcr  the  latter  stages  of  loading  with  more  dense  moire  fringe  patterns. 

Figure  6  shows  the  various  J-values  with  crack  extension  computed  by  the  build  in 
algorithm  of  the  commercial  FEM  code  used,  Jfem,  and  through  contour  integration  [6,7,8] 
using  the  moire  fringe  data,  Jexpl  and  Jexp2-  Jfem  computation  had  to  be  terminated  at  Aa  = 

1 .5  mm  due  to  the  large  distortion  of  the  small  finite  element  surrounding  the  crack  tip.  The 
equivalent  elastic  J-value  was  obtained  by  setting  J  =  G  (strain  energy  release  rate),  Jlefm. 
and  Jshih  was  computed  using  the  procedure  and  tables  in  [9,10].  Jexpl  and  Jexp2  differ  by 
the  size  of  contour  with  former  being  5x5  mm  and  the  latter  was  10  x  10  mm.  The 
differences  in  the  theoretical,  experimental  and  numerical  results  essentially  follow  the 
same  trend  in  [6,8J.  It  is  interesting  to  note  that  both  JpEM  and  Jshih  are  relatively  close 
and  both  differ  substantially  with  Jexpl  and  Jexp2- 

Figure  7a  shows  log-log  plots  of  the  measured  v-displacements,  the  computed  HRR 
and  the  I Q I  components  of  the  v  displacements  along  the  radial  line  of  0  =  45°  for  a  stable 
crack  growth  of  Aa  =  1.5  mm.  Here  I Q I  is  the  absolute  value  of  the  difference  between  the 
experimental  and  HRR  displacement  components  and  is  not  the  Q  value  in  Equation  (5). 
The  measured  v  displacements  and  the  HRR  components  of  the  v  displacements  nearly 
coincide  for  rOo/J  >  10  but  differs  substantially  for  rOo/J  =  1-  The  I Q I  component  of  the  v 
displacements  is  an  order  of  magnitude  smaller  and  possibly  within  the  error  band  of  the 
measurement. 

Figure  7b  shows  the  corresponding  log-log  plot  of  the  measured  u-displacement  and 
the  computed  HRR  and  the  I Q I  components  of  the  u  displacement.  For  rOo/J  =  10  the  I  Q ! 
components  is  the  same  order  of  magnitude  as  the  measured  u  displacement  and  the 
corresponding  HRR  component  and  varies  approximately  as  rl  l874  for  r  =  l  mm.  This 
exponent  indicates  that  the  corresponding  normal  strain  component  is  nearly  a  constant 
independent  of  r  and  is  closer  to  the  T-component  as  predicted  by  the  J-T  theory  [20], 
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Figures  8a  and  8b  show  the  log-log  plots  of  the  measured,  computed  HRR  and  Q 
component  of  the  v-  and  u-displacements,  respectively  for  a  large  crack  extension  of  Aa  = 
3.9  mm.  While  the  HRR  and  measured  v-displacements  for  njo/J  =  10  are  in  reasonable 
agreement  despite  the  large  crack  extension,  significant  difference  is  noted  in  the  u- 
displacement.  The  associated  Q  component  is  negligible  for  the  v-displacement  but 
increases  by  1^-34  for  the  u-displacement.  Thus,  for  this  larger  Aa  =  3.9  mm,  neither  the  |-Q 
nor  the  |-T  theories  appears  to  be  applicable. 

DISCUSSION 


The  uO),  e(l)  and  functions  are  yet  to  be  determined  for  the  state  of  plane  stress 
assuming  that  the  stresses  are  separable  into  functions  of  r  and  0.  The  negligible  !  Q  I  - 
component  of  the  v-displacement  and  the  slope  of  the  I Q I  -component  of  the  u- 
displacement  for  small  stable  crack  growth  suggest  that  the  J-T  theory  may  account  for  the 
cummulated  effects  of  the  higher  order  terms  in  Equation  (2). 

CONCLUSIONS 

A  comparison  of  experimentally  determined  crack  tip  displacement  fields  in  a  thin 
SEN  specimen  with  predictions  by  the  J-Q  theory  has  been  carried  out.  The  comparison  is 
qualitative  since  the  form  of  the  J-Q  displacement  fields  is  not  known  for  plane  stress. 

Limited  experimental  results  involving  the  crack  tip  displacement  fields  in  thin 
aluminum  SEN  with  small  stable  crack  growth  showed  that  the  J-Q,  based  on  the  plane 
strain  form,  may  not  be  present.  On  the  other  hand,  the  simpler  J-T  crack  tip  field  [20] 
could  exist. 

Also,  the  use  of  the  present  J-integral  computation  procedures  using  [9]  and  a 
commercial  code  for  elastic-plastic  fracture  analysis  of  thin  plates  in  the  presence  of 
stable  crack  growth  must  be  reinvestigated  in  view  of  the  observed  large  discrepancies. 
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Fig.  2.  SHN  2024-T3  aluminum  alloy  specimen  with  J-inlegral  contours. 
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Fig.  3  Finite  element  mesh  used  in  elastic-plastic  analysis. 
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Fig.  4.  Applied  load  vs.  load  line  displacement  for 
2024-T3  SEN  specimen. 
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Fig.  6.  J-integral  for  2024-T3  SEN  specimen. 
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Fig.  7a.  v-displacemenl  in  2024  T3  SF'N  specimen. 
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u-displacemont  in  2024-13  SliN  specimen. 
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